Entanglement and boundary theories of Resonating Valence Bond wavefunctions 
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We investigate the relations between bulk and boundary properties of short-range resonating 
valence bond wave functions using Projected Entangled Pair States. We show that, for a bipar- 
tition, the boundary Hamiltonian defined on the edge can be written as a product of a highly 
non-local projector, which fundamentally depends upon boundary conditions, with an emergent 
(local) one-dimensional (superfluid) t-J Hamiltonian, which arises due to the symmetry properties 
of the auxiliary spins at the edge. This multiplicative structure, a consequence of the disconnected 
topological sectors in the space of dimer lattice coverings, is then characteristic of the topological 
nature of the states. For systems with open boundary conditions, it is shown that the entangle- 
ment spectrum, which reflects the properties of the (gapped or gapless) edge modes, is a subset of 
the spectrum of the local Hamiltonian, also explaining the origin of the topological entanglement 
entropy. We propose to use these features to probe topological phases in microscopic Hamiltonians. 
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Conventional ordering in condensed matter systems 
is traditionally associated to symmetry-breaking and to 
the existence of a local order parameter (Landau the- 
ory). Topologically ordered phases of matter [1 offer 
completely new classes of systems for which the ground 
state (GS) degeneracy depends on topology (disc, cylin- 
der, torus, etc.). The (short-range singlet) Resonating 
Valence Bond (RVB) wavefunction proposed by Ander- 
son [2] as the parent Mott insulator of high-temperature 
superconductors is a celebrated example. Such topologi- 
cal phases carry emerging fractionalized excitations and 
raise growing attention due to their potential to realize 
fault-tolerant setups for quantum computing. [3 

Interest for topological liquids in quantum antiferro- 
magnets [J and in related microscopic models [5-8 has 
also triggered the search for novel theoretical tools ca- 
pable of detecting topological order, in particular en- 
tanglement measures used in quantum information. A 
common setup consists of dividing the system into two 
regions (named A and B) and compute the reduced den- 
sity matrix (RDM) in the GS of e.g. the A subsystem. 
The entanglement entropy (EE), defined as the Von Neu- 
mann entropy of the RDM S'vn = —pA In p^, contains an 
extensive term - proportional to the length of the bound- 
ary (area law) - and a universal sub- leading constant, the 
topological EE. Specific disc-like setups [9 or cylindrical 
geometries can be used to extract the topological EE. 

In fact, — In pA can be seen as a (dimensionless) Hamil- 
tonian i^5, a key conceptual object. First, its spectrum, 
the so-called entanglement spectrum (ES), has been con- 
jectured to show a one-to-one correspondence with the 
spectrum of edge states. This remarkable property was 
first established in fractional quantum Hall states [10] 
and, then, in quantum spin systems |lT]. Furthermore, 
Projected Engangled Pair States (PEPS) flT offer a nat- 



ural formulation of the relation between bulk and bound- 
ary. In [13 , an explicit isometry was constructed which 
maps the Hamiltonian H^j onto another one acting 
on the space of auxiliary spins living at the edge of re- 
gion A, while keeping the spectrum. Furthermore, for 
various two-dimensional (2D) models displaying quan- 
tum phase transitions, like a deformed AKLT [14] or an 
Ising-type [15 model, it was found flT that a gapped 
bulk phase with local order corresponds to a boundary 
Hamiltonian with local interactions, whereas critical be- 
havior in the bulk is reflected on a diverging interaction 
length of Hi). 

Entanglement properties of 2D topological phases are 
less well understood. Rokhsar-Kivelson (RK) wave 
functions, defined as equal- weight superposition of fully 
packed dimer coverings, exhibit critical behavior on bi- 
partite lattices [16 or realize the simplest topologi- 
cal phase, the so-called Z2 liquid, on frustrated lat- 
tices. [17 The topological EE of critical and topologi- 
cal RK wavefunctions have been computed using vari- 
ous topologies [18 and the boundary Hamiltonian corre- 
sponding to the GS of Kitaev's toric code [3 was shown 
to be non-local. [13] Unfortunately, RK-like wavefunc- 
tions are not generic - their ES is completely dispersion- 
less [19] - and do not describe real quantum S = 1/2 
spin systems. [20l [21] In contrast, (short-range) RVB 
states, defined as linear superposition of hardcore cov- 
erings of non- orthogonal nearest-neighbor SU(2) singlets 
(see Fig. [ija)), appear to be closer to physical sys- 
tems. Very recently, the (Renyi) EE between (finite) 
cylindrical regions has been computed [22 numerically 
for the critical [23l [25 RVB state on the square lat- 
tice. Similarly, the (Renyi) topological EE of SU(2)- 
symmetric gapped chiral and Z2 spin liquids was ob- 
tained [25 using Kitaev-Preskill prescription. Neverthe- 
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FIG. 1: (Color online) (a) Particular dimer hardcore covering 
on a Nh X Ny cylinder. Ellipses represent singlets of two spins 
1/2. The RVB wavefunction is defined as the equal- weight 
superposition of all such configurations. Bi-partitioning the 
cylinder defines two (virtual) "L" and "R" edges, (b) PEPS 
on the same cylinder involving identical (rank-5) tensors (blue 
squares) on the sites with four virtual indices {D — 3) and one 
physical variable (curly line). The "Bl" and "Bi^" boundaries 
are defined by the choice of the external indices (OBC corre- 
spond to choose "2" for all). Here, we take the limit Nh oo 
(infinitely long cylinder). 



less, ES and boundary Hamiltonians of such RVB/spin 
liquids wavefunctions is unknown. In this Letter, we 
study entanglement properties of both critical (square 
lattice) and gapped topologically-ordered (kagome lat- 
tice) RVB wavefunctions [24 on infinite cylinders mak- 
ing use of simple PEPS representations. In particular, we 
show that the boundary Hamiltonian can be written as 
a product of a highly non-local projector, which depends 
fundamentally on the boundary conditions, by a local 
one-dimensional t-J model, which arises due to the sym- 
metry properties of the auxiliary spins at the boundary 
and characterizes the (gapless or gapped) edge modes. 
This multiplicative structure, a consequence of the dis- 
connected topological sectors in the space of dimer cov- 
erings of the lattice, reflects then the topological nature 
of the states. 

PEPS representation of RVB states. - We start with 
the square lattice RVB wavefunction (NN | ti) - I It) 
singlets are all oriented from one sublattice to the other) 
on a cylinder of length and circumference N^^ de- 
picted in Fig.jlja), which can be expanded in the local Sz- 
basis, I^rvb) = c^l^i, S2, . . . , sm), where Sn = 0, 1 
are qubits (representing the two Sz = ±1/2 spin com- 
ponents) on the M = N^Ny sites and S = {sn}- Such 
a state can be represented by a I) = 3 PEPS [151 EI] 
(up to local unitaries) where each lattice site is replaced 
by a rank-5 tensor A^a^a';(3,(3' labeled by one physical in- 
dex, s = or 1, and by four virtual bond indices (vary- 
ing from to 2) along the horizontal (a, a') and verti- 
cal (/3, /3') directions. Physically, the absence of singlet 
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FIG. 2: (Color online) (a) Kagome lattice. The L and R edges 
are shown (thick dotted line). Three site tensors (red dots) 
carrying the physical indices and two 120-degree tensors (in 
the center of the shaded triangles) are grouped together (a,b) 
to construct the basic tensor (c). 



on a bond is encoded by the virtual index being "2" on 
that bond. To enforce the hardcore dimer constraint, 
one takes ^/ = 1 whenever three virtual indices 

equal 2 and the fourth one equals s, and /3' = 

otherwise. The amplitudes cs are then obtained by con- 
tracting all virtual indices, except the ones at the ends 
of the cylinder fixed by boundary conditions (BC), as 
depicted in Fig. [TJb). For the kagome lattice, as shown 
in Fig. [2j the RVB state can be represented in terms of 
rank-3 tensors, (i) A^.^ on the sites - A^.^ = = 1 
and zero otherwise - and (ii) on the center of each trian- 
gle, i^2,2,2 = 1, and Ra,i3,j = ^ai3-f Otherwise, with eaf3j 
the antisymmetric tensor. [21] One can then group the 3 
sites on each unit cell to obtain a rank-5 tensor (the phys- 
ical dimension is now 2^ = 8) connected on an effective 
square lattice (Fig. [2|b,c)). Note that for the Kagome 
PEPS, one can find a local parent Hamiltonian for which 
the degeneracy is equal to 4 on the torus. [21] 

RDM on infinite cylinders - The boundary Hamilto- 
nian Hi) can be derived from the reduced density operator 
cr^ = exp {—Hb) acting on the edge indices, following the 
procedure given in [13 . For the kagome lattice, there is 
no reflection symmetry w.r.t. the cut so the RDM for 
the left (right) side takes the form = \/^crL\/^ 
(cr^^ = y^o-Ry^) where (Jl and are obtained by 
contracting the tensors of the left and right half-cylinders, 
respectively. [13 Note that <j^^ and <j^^ give identical 
ES. Here we are interested in RVB cylinders with infi- 
nite lengths in both directions. First, we fix the cylinder 
perimeter {Ny = 4, 6, 8) and take the limit A^^ — ?► oo (in 
practice, the RDM for = 50 is fully converged). 

Disconnected sectors - A crucial feature of topologi- 
cal states is that the RDM depends intrinsically on the 
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FIG. 3: (Color online) ES (w.r.t. the GS energy ^o) of an 
infinitely long kagome cylinder of perimeter Nv — 8. Eigen- 
states with half-integer (a) and integer (b) spins correspond 
to odd and even sectors, respectively (see text), (c) Weights 
\dx^{r)\'^ of the 2-body operators appearing in i/iocai- 

choice of the BC {Bl = for simplicity), even when 
Nh oo. Indeed, due to local constraints, (vertical) con- 
figurations of virtual indices can be grouped in disjoined 
conserved sectors which can be addressed independently 
from proper choices of Bl = Br. Hence, although the 
RDM (and Hb) acts on all 3^^ degrees of freedom of the 
L (or R) edge, in each sector it contains a finite frac- 
tion of zero-weight eigenvalues, i.e. a finite fraction of 
eigenstates of ^5 have infinite-energy. Calling V the 
projector on the finite-energy subspace, we split Hjj as 
Hb = Hi^ Aopo(Id - V). where Id is the 3^- x 3^- 
identity-operator, Ptopo ^ co, and Hi is supported by 
the non-zero eigenvalues sector of the RDM. More pre- 
cisely. Hi can be factorized as Hi = H\oca\V where H\oca\ 
is a Hamiltonian (shown later to be local) acting on the 
whole boundary space [28 , commuting with V and in- 
dependent on BC. The Kagome cylinder has only two 
sectors defined by 7^ = Veven and V = Vodd = Id — Peven , 
which can be obtained by choosing an even or odd num- 
ber of "2" external (virtual) indices for Bl = Br, respec- 
tively. For the RVB wavefunction on the square lattice, 
there are Ny -\- 1 disconnected sectors defined by projec- 
tors Va enforcing a fixed difference A^2,x — A/'2,y = A of 
the numbers of "2" on two X and Y alternating sublat- 
tices on the edge, A = —Ny/2, • • • , Ny/2. The fact that 
Hi is known for all sectors implies that i^iocai is uniquely 
determined. [29] Open BC (OBC), relevant for physical 
systems, are defined by having "2" set on all the Bl and 
Br indices (Fig. [ijb)). Therefore OBC ffive access to 
the even (A = 0) sector on the kagome (square lattice) 
cylinder. [30^ 

Bipartite entanglement spectrum - The full ES is given 



FIG. 4: (Color onUne) Same as Fig. [s] for the critical RVB 
state (square lattice). Eigenstates with half-integer (integer) 
spins correspond to A odd (even) - see text. For OBC one 
gets a subset of (b) (A = sector), (c) Weights \dxfx{r)\'^ 
of the 2-body operators appearing in ifiocai- The diagonal 
interaction ^2X2^^ shows a long-range behavior. 

by the spectrum of i^iocai- Our results are summarized 
in Figs. |3ja,b) and [4]^a,b), for infinitely- long cylinders 
with kagome and square lattices. The (excitation) ES are 
shown as a function of momentum around the cylinder 
and the eigenstates are labelled according to their spin- 
multiplet structure inherited from the su(2) symmetry 
of the RVB state, although with the 1/2 representa- 
tion. [26 A careful finite size scaling suggests that the 
kagome (square) lattice cylinder ES is gapless (gapped) 
in the limit Ny ^ 00. Since these features are opposite 
to what is expected in the corresponding bulk systems, 
we deduce that the ES characterizes specifically the na- 
ture of the L and R edge modes (Fig. [T]). Note that for 
given choice of BC, the actual ES is the spectrum of a 
projected Hi Hamiltonian and, hence, is a subset of the 
full ES. 

Boundary Hamiltonian and non-locality - Next, we 
investigate the exact connection between the boundary 
Hamiltonian and the bulk and topological properties of 
the system. We look for its explicit form, trying to make 
the connection with D = 3 models, with su(2) symmetry 
corresponding to the 1/200 representation. Any opera- 
tor Oedge acting on the edge can be expanded in terms of 
^2N^ orthogonal operators. For this purpose, we use a lo- 
cal basis of 9 (normalized) operators {xq, • • • ,^8} which 
act on the local basis of configuration (at some site z), 

{|0),|1),|2)}, e.g. xo = 1, ii = vf(|0)(0| - |1)(1|) 
and X2 = 73(|0)(0| + |l)(l| - 2|2)(2|), for the diag- 
onal matrices, complemented by X3 = x\ = -\/3|0)(l| 
acting as "spin" operators, and x^ = xl = -\/3|2)(0| 



and Xq 



_ ;9.t _ 



>/3|2)(l| acting as annihilation and cre- 
ation (hardcore) bosonic operators. [27| The expansion 
in terms of N-body operators reads, 



Oedge = CoNy 



E 



eA^.(r,r')4<+'-xt+'^'+---, (1) 



where each group of terms involves products of (1 < 
N < Ny) on-site xx (A 7^ 0) operators. Here the sums are 
restricted to non-equivalent relative distances and only 
translations giving distinct sets of sites are performed. As 
seen from the distribution of their weights in Fig. |5|a), 
projectors are highly non-local, conferring a fundamen- 
tally non-local character to the boundary Hamiltonian 
Hi. This is also to be expected for realistic topological 
GS of microscopic Hamiltonians on geometries involv- 
ing open or fixed BC in some directions. In contrast, 
^locai [29" contains primarily one- and two-body contri- 
butions (Fig. (sja)). The weights of all (non-zero) two- 
body terms of i^iocai shown in Figs. [3]^c)j4jc) reveal ex- 
ponential decays with distance r, except for a diagonal 
term in the case of the square lattice. These remarkable 
features are to be connected to the bulk correlations of 
the RVB wavefunctions, short-range in the kagome lattice 
and critical for the square lattice. Interestingly, i^iocai is 
formally equivalent to a one-dimensional "t-J model" 
describing motion of "holes" (the "2") in a spin fluctu- 
ating background (the '0" and "1" qbits) supplemented 
by additional density-density and pair-field terms, con- 
ferring a superfluid character to the edge. 

Topological entropy. - For infinite cylinders of the 
kagome lattice with OBC, the EE Svn = — Tr{cr^ Ina^} 
can be well fitted according to Svn = Sq -\- ANy^ as 
shown in Fig. jsjb), where 5*0 is the topological EE. 
The existence of a finite 5*0 is a direct consequence of 
the particular structure of Hi. Indeed, the EE is given 
(crudely) by — In A/* where A/" is the number of eigenstates 
of Hi below a fixed energy scale of order 1. For OBC, 
the support of Hi = i^iocai^^even coutaius A/" — ^ 3^^ 
states and = — In 2, as expected for a topological Z2 
spin liquid. Summing over all sectors amounts to take 
Hi = i^iocai so that all eigenstates of the ES contribute 
and Svn oc Ny (as can be shown rigorously) as seen in 
Fig. jsjb). For the square lattice, severe constraints leads 
to an extensive number (i.e. proportional to the perime- 
ter Ny) of disconnected sectors on the edge of dimension 
A/" ^3^^, therefore introducing negative logarithmic 
corrections ~ —InNy to the EE for any BC (see e.g. data 
for OBC on Fig. [sjb)). The long-range diagonal interac- 
tion in i^iocai (Fig.[4]^c)) is also responsible for deviations 
from the area law, even when considering all sectors. 

Discussion and outlook. - We have shown that bound- 
ary Hamiltonians associated to gapped and gapless RVB 
states can be written as i^iocai7^ + Aopo(Id — 7^), /3topo ^ 
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FIG. 5: (Color online) (a) Distribution of weights vs N (semi- 
log plot) of T^even, Vo " coiiesponding to OBC - and of i/iocai 
for kagome (circles) and square (squares) lattices (Ny — 8). 
(b) Entanglement entropy versus Ny for OBC (open symbols) 
or when summing over all sectors (shaded symbols). 



00. In particular, we have established the existence of 
a projector V (which intrinsically depends on the BC) 
onto a restricted subspace at the edge (as for Kitaev toric 
code [3 , 13 ), a consequence of the disconnected topolog- 
ical sectors in the space of dimer coverings of the lat- 
tice. We argue that the non-local character of the result- 
ing boundary Hamiltonian is the fingerprint of topology. 
The ES is a subset (associated to V) of the spectrum of 
the emerging local Hamiltonian i^iocai acting on the un- 
restricted edge space. In contrast to the toric code for 
which i^iocai is trivial, here i^iocai takes the form of a 
short-range (bosonic) t-J model (including a long-range 
diagonal interaction for the critical RVB state). We pro- 
pose to use the features of the boundary Hamiltonian to 
detect topological phases in microscopic models. [6^-^ 
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